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Abstract: In nonlinear physics, the interactions among solitons are well studied thanks to the multiple 
soliton solutions can be obtained by various effective methods. However, it is very difficult to study 
interactions among different types of nonlinear waves such as the solitons (or solitary waves) , the cnoidal 
periodic waves and Painleve waves. In this paper, the nonlocal symmetries related to the Darboux trans- 
formations (DT) of the Kadomtsev-Petviashvili (KP) equation is localized after imbedding the original 
system to an enlarged one. Then the DT is used to find the corresponding group invariant solutions 
such that interaction solutions among different types of nonlinear waves can be found. It is shown that 
starting from a Boussinesq wave or a KdV-type wave, which are two basic reductions of the KP equation, 
' the essential and unique role of the DT is to add an additional soliton. 
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1 Introduction 

To study the nonlinear waves, there are typical important prototype models such as the KdV (Ko- 
rteweg de-Vries), KP (Kadomtsev-Petviashvili), NLS (nonlinear Schrodinger) and SG (sine-Gordon) mod- 
els. These systems appear in almost all the physics fields, especially, in fluid physics, plasma physics, 
optics, condense matter, quantum fields and astrophysics pQ. Many kinds of nonlinear waves (nonlinear 
excitations) such as the solitons, conoidal periodic waves, Painleve waves (solutions of the Painleve tran- 
scendent I- VI) are found by using various effective methods. However, to study the interactions among 
these nonlinear waves is still a thorny issue though the interactions among solitons have been studied 
well. 

There are some excellent approaches to find multiple soliton solutions of integrable system. For 
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instance, the simplest method may be the Hirota's bilinear method [2]. Usually, whence a nonlinear 
system is changed to a bilinear form, its multiple soliton solutions can be directly written down and then 
the interactions among solitons can be studied in detail from the analytic expressions. The Darboux 
transformations (DT) and/or the Backlund transformations (BT) are two other effective methods [3] . 
Using these two methods, in principle, one can obtain a new solution from a known one. However, in 
practice, one can only find multiple soliton solutions starting from simple constant solutions. It is still 
very difficult to find new explicit solutions starting from nonconstant nonlinear waves such as the cnoidal 
waves and Painleve waves via DTs and BTs. It is not very clear that what kind of solutions can be find 
via DT and/or BT if one takes a nonconstant and non-soliton solutions as a seed. 

Since the Lie group theory was introduced by Sophus Lie to study differential equations [4] , the study 
of Lie group has always been an important subject in mathematics and physics. Using both classical 
and non-classical Lie group approaches 16] , one can reduce dimensions of partial differential equations 
(PDEs) and proceed to construct the analytical solutions of these PDEs. 

Traditionally, to study symmetry groups one studies their symmetry algebras first and then the 
related groups via Lie's first theorem because the algebra problems are linear while the group problems 
are nonlinear. The recent studies show us that for many kinds of nonlinear partial differential equation 
systems, to firstly find symmetry groups of the nonlinear systems may be much more convenient than to 
firstly find related symmetry algebras [7]. 

Actually, there exist some other types of examples in soliton theory in which the groups (finite trans- 
formations) such as the DTs and BTs have been successfully studied while the corresponding symmetry 
algebras are not very familiar. In Refs. [H [9], the authors have shown us that the symmetries related to 
DTs for some well known equations, such as the KdV equation, the KP equation, the CDGKS equation, 
etc. are nonlocal. It is known that using Lie point symmetries, the dimensions of the related PDEs can 
be reduced. An important problem is how to use the DT related symmetries to reduce the dimensions 
of the original PDEs? A simple direct way to solve this problem is to prolong the system such that 
the nonlocal symmetry of the original system becomes a local one of the prolonged system |10j . In this 
paper, we will show that whence the DT related nonlocal symmetries are localized, it will be quite easy 
to obtain new exact solutions from non-soliton type seed solutions via DT related symmetry reduction 
approaches. 

In Sec. 2 of this paper, the first DT related nonlocal symmetry of the KP equation is localized. In Sec. 
3, using the localized symmetry of the first DT, the Levi transformation (the second DT) is re-obtained via 
Lie's first theorem. In Sec. 4, the similarity reductions for the prolonged system are considered according 
to the standard Lie point symmetry approach. Some corresponding explicit interaction solutions among 
solitons and Boussinesq waves and KdV type waves are given in Sec. 5. The last section is devoted to a 
short summary and discussion. 
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2 Localization of the non-local symmetry related to the first DT 
of the KP equation 

We consider the KP equation [TT] in the form 

(u t - 6uu x + u xxx ) x + 3u yy = 0, (1) 

where the subscripts denote derivatives. The KP equation was firstly derived to study the evolution of 
long ion-acoustic waves of small amplitude propagating in plasmas under the effect of long transverse 
perturbations [TT]. The KP equation was widely accepted as a natural extension of the classical KdV 
equation to two spatial dimensions, and was later derived as a model for surface and internal water waves 
|12j . and in nonlinear optics |13] and almost in all other physical fields such as in shallow water waves, 
ion-acoustic waves in plasmas, ferromagnetics, Bose-Einstein condensation and string theory. The KP 
equation is also used as a classical model for developing and testing of new mathematical techniques. 
It is well known that the KP equation possesses the Lax pair 

^XX - + tPy = 0, (2) 

+ H>xxx - 6uijj x - 3 I u x — i Uydx ) if) = 0. (3) 



and the adjoint Lax pair Eq.([T]) 

4>XX - Ucj) - <j>y = 0, (4) 



4>t + ^<(>xxx - Gu<f> x - 3 yu x + j Uydxj = 0. (5) 

That is to say the integrable conditions of Eqs. (l2l3p and (|4l5p . ip y t = i/jt y and 4>yt = 4>ty, are just the KP 
equation dJJ). 

A symmetry, a u of the KP equation is defined as a solution of its linearized equation 

(af - 6a u u x - dual + a u xxx ) x + 3a^ = (6) 
that means ([TJ is form invariant under the transformation 

u + ea u (7) 

with the infinitesimal parameter e. 

In Ref. [8], the non-local symmetry corresponding to the DT of the KP equation [14] has been 
derived, which is in the form 

° u = WOx, (8) 



where ip, <ft satisfy the Lax pairs (121.'] I) and (|4I5I) . 

Because (O is just the infinitesimal form of the DT, one can naturally believe that it can be localized 
to a Lie point symmetry such that we can use the Lie's first theorem to recover the original DT. From 
|8"]). it can be apparently seen that the non-local symmetry contains the space derivative of functions 
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ip and (j). Then, to localize the nonlocal symmetry (|8|), we have to know what are the corresponding 
transformations for the quantities ip, <$>, 

V>i=V>», (9) 

and 

4>i = 0, (io) 

whence the field u has the symmetry transformation u — > u + ea w . In other words, we have to solve the 
linearized equations of © , ©, © and (p~0|) . 

at - a^ - uo* + a* = 0, 

a* = a* (11) 

whence a u is given by ([5]) . 

It is not difficult to verified that the solution of (jTTJ) with © has the form 

^ = Tj&P, 
1 



a 

a* 1 = \{^iP + 4> 2 ^), (12) 

where the new quantity p is defined as 

Px = ip<t>- (13) 
Two compatibility conditions of (|13l) are worth to mentioned here 

P v = "001 - ipt<f>, (14) 
p t = 4^i0i + 6uijxt> - 4(fyipix - 4ipcj)i x , (15) 

that means the conditions p x t = ptx , Pty — Pyt and p xy = p yx are satisfied identically. 

The appearance of the quantity p in the symmetry solution (|12[) . we have to further solve the 
linearized equation of (|13p 

a P = a^tp + V-cr* (16) 

with the condition (|12[) . 

It is interesting that the solution of (|16p with (jT2"j) has the simple form 

aP = \p\ (17) 
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The result (fl7|) hints us that p is a solution of the Schwarzian KP equation, 



-I X 



S(gj] =0. (IM 



where the Schwarzian derivative {p; xj = — S^ay and the quantities — and — are all invariant 
under the Mobious transformation invariant with the infinitesimal transformation (|17|) . This fact implies 
a new method to find Schwarzian forms of nonlinear systems which may be more significant in the future 
studies in discrete systems where the traditional Paileve analysis [15] and the singularity manifold do not 
work. 

The results (fT2j) and ([TT]) show us that the DT related nonlocal symmetry ([8]) in the original space 
{x 7 y, t, u} has been successfully localized to a Lie point symmetry in the enlarged space {x, y, i, u, ip, cj>, 
ipi, (f>i, p} with the vector form 

V = (V^i + Hi)d u + i^pfy + ^cf>pdt + i(^ lP + ^^d.^ + hfap + 2 V)<9 01 + ^p 2 d p . (19) 

In other words the DT related symmetry is only a special Lie point symmetry (|19p of the equation system 
{©, ©, CP, ©, (HOj), (H3J)}. The equations ©, ©, (HU) and ([15]) may should also be included in the 
enlarged system. However, we only treat them as the compatibility conditions of (HK2), (Hill]) and so on. 



3 Second type DT of the KP equation 

In this section, we study the finite transformation form of (|19[) by using Lie's first theorem. To 
find the nonlocal symmetry ([8]) from DT, we take the infinitesimal group parameter is just the spectral 
parameter [Sj |S]. So it is natural that the finite transformation of (fTO)) is just the usual first type of 
DT if one take the spectral parameter as group parameter. Here we are more interest to ask if we take 
the group parameter is independent of the spectral parameter, what kind of finite transformation will be 
found? That means we should solve the following "initial value problem" 

du'(e) 



de 



= (^i + <Hi) u'(0) = 



de 2 

^ = \<t>'{e)p>{e), 0'(O) = 0, 

= i W(e y (e) + V) ' (e) V'(e), 1&J(0) = ^, 

d#(e) 1 



de 2 



(^(cy(e)+^(e)V(e), #(0) = 



^ = ^W a . *i(°)=i». (20) 

with the group parameter e being independent of the spectral parameter. 

By solving the initial value problem (|20[) . it is not difficult to find the following theorem. 
Theorem 1. If {u, ip, </>, ipi, <fii, p] is a solution of the extended system ([l])-©, ((9J, (fT()|) and (fT3|) . 
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then so is {«'(£), 0'(e), 0'(e), ^(e), <f>[{e), p'(e)} with 

2p 



2 - ep' 

2^ 
2 - ep' 

20 
2 - ep' 
2V>i 2c0 2 



2 - ep (2 - ep) 2 ' 

20i 2e0 2 '0 
2-ep + (2 - ep) 2 : 



2e(V>0i + 0^i) 2eW 
U(£) = " + 2-ep + (2~^ (21) 

for arbitrary group parameter e. 

It should be emphasized that this finite transformation theorem is not the usual DT (the first type 
of DT). It is equivalent of the so-called Levi transformation, the second type of Darboux transformations 
|16j . In other words, at the algebra level, two types of DTs share a same form but with different group 
parameters, the first one is just the spectral parameter while the second one is an independent parameter. 

From the first equation of ([2"Tj) . we know that the second DT is nothing but the Mobious transfor- 
mation of the Schwarzian forms of the original nonlinear systems. 



4 New symmetry reductions of the KP equation 

Now we take into account the Lie symmetries of the closed enlarged system. It is known from the 
standard Lie point symmetry group approach that, in order to find the Lie point symmetry, we may 
assume the symmetries have the vector form 

d d d d d d d d d , 

V = X \-Y \-T \-U h* h$ h*i h$i VP — , (22a) 

dx dy dt du dtp 90 dtpi 90i dp 

where X, Y, T, U, "f, $, $i, P are the functions with respect to {x, y, t, u, ip, 0, ipi, 0i,p}, which means 
that the closed system is invariant under the transformations 

{x,y,i,u,^,0,^l,0i,p} -4 {a; + eX,y + ey,t + eT,7i + e[/,^ + e*,0 + e$,Vi + e*i,0i+e$i,P + e-P} (23) 

with a small parameter e. Equivalently, due to the fact that the system is not explicitly space-time 
dependent, the symmetries in the vector form (|22|) can be written as a function form 



<j u = Xu x + Yu y + Tu t - U, (24a) 

^ = Xip x + Yip y + Tip t - (24b) 

(T+ = X<j> x + Y<}> y + T0 t - (24c) 

a in = X4> lx + Yi^ lv + Tipu - *i, (24d) 

ofr = X(f> lx + Y<f> ly + T<p lt - $i , (24b) 
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a p = Xp x + Y Py + Tp t - P. (24f) 

In this notation, er u , tr^,, a^, a i p 1 , a^ >l , a p are the solutions of the symmetry equations, i.e., the linearized 
equations for the closed system 

&u,xt X2iU x (7 u Qu(T u xx Q&uUxx ~t~ @u,xxxx ~t~ 3(T u ^yy 0, (25a) 

cr^,xx - ua^, - <7 u ijj + a^, : y = 0, (25b) 

<J<p,xx - ua^ - a u (j) - a^ y y = 0, (25c) 

cr^ 11 - <r^,x = 0, (25d) 

0>i - (T^, x = 0, (25e) 

o-p - ipa$ - a^<f> = 0. (25f) 

Substituting Eq. (|24[) into Eq. (l25[) and eliminating u yy , ip X x, 4>xxi i>x, 4>x and p x in terms of the closed 
system, we get more than 500 determining equations for the functions X, Y, T, U, ^f, $, $i and P. 
Calculated by computer algebra, the general solutions of them take the form 

X = \f t x - ^-fttV 2 -~thV - 9, Y=\f t y + h 1 T = f, 
3 18 6 3 

2 1 1,1 1 

U = -^ftu-—fttx + —fttty+^huy + -gt-c^i(f>-c(pi-ip, 

3 18 108 36 6 

* = {^ fttty3 -h xftty -^ xht + ^ htty2 + l 9tV+F 'Y)^ 

* = (j^xfttV + J^ xh t - T^fttty 3 - -f^ttV 2 ~^9tV - F - i/ t + c 2 - y 
*i = (^^futy 3 - ^xh t ~ ^xfuy + ^h u y 2 + ^g t y - ^ft + F - -tp! 

Y8 ftty + T2 ht + C -¥)^ 
$1 = (~^^**y 3 + ^- t/i * + j$ x ftty - ^ h tty 2 - \gty - |/t - ^ + c2 - y 

J 3 = {\ft + o 2 )p + c 1 - C -^- 1 (26) 
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where {/, h, g, F} are arbitrary functions of t and {c, ci, C2} are arbitrary constants. Consequently, it is 
convenient to rewrite symmetries (1241) as 

1 1 2 1 \ (2 \ 2 1 1 2 

■^ftx - —ftty - -h t y — g\ u x +\ -ftV + hju v + fu t + -f t u + j^fttx - —fmy 

■^ h tty - \gt + c(#i 
36 6 

( ~ft% - j^fuy 2 - ^fky -g^j4> x + (\ftv + hjip y + fip t + {r^ x fttv + ~& xht ~ ^fttty 3 
^2 h tty 2 - ^gty - Fj ip + ^cpip, 

\ftx - -^ftty 2 - ^hy -g^(j) x + \\ftV + M fa + ffa ~ \T$ x ftty + J2~ xht ~ ^ftttV 3 

1 2 1 1 1 

■^ h tty - ^gty - F - -f t + c 2 - -cp 

>>■, = \\ftx - j^ftty 2 ~ ^hy - g^j fa x + f|/t2/ + hj tfj ly + ftpu + (J^ftty + J^ h ^j i> 

f 1 . , 1 , 1 . 1 , o 1 1 . „\ , 1,1 



1 32 ^fttty - yi xh t ~ jgxftty + j2 h tty + ^gty - g/t + F I + -^cp^x + ' 

g/tz - jgftty 2 - ^h t y - g^j fa x + f ~/t2/ + hj fa y + ffa t - (jgfttV + j^ h t 

^fttty 3 - j2 xh t ~ jgZftty + Y2~ h tty 2 + ^gty + ^ft + F- c^j fa + ^cpfa + ^cfa^j, 

- j^ftty 2 ~ ^h t y -g^Jp x + (\ftV + h) p y + fp t - c 2 p - c x + ^cp 2 . (27) 

On observation of (|2~T1) , it is not difficult to find that the c-dependent generator of the related Lie symmetry 
group of Eq.([T]), (facfn + faip)-§-, is just that corresponding to the Darboux transformation. 

Similarity solutions to the closed system can be found by setting ct m , cr,/,, er^, er^ , cr^ , <7 p be zeros in 
Eq. (|27l) . which is equivalent to solving the characteristic equations 

dx dt du dip d(j) dfa dfa dp 

In c = case, three reduction cases have been studied by many authors, say, in [17 . In the proceeding 
paragraphs, three new cases which are generalizations of the known three reductions for c ^ are 
discussed in detail. 

Case 1. / ^ 0. In this case, without loss of generality, we rewrite the arbitrary functions {g, h, F} in 
Eq.||2I|) as 

h = yotfK h= fr ^t, 



9 = qJ Vot -J ' x ot, 
F = ^y 3 otf 3 + fitf + \yof(x tf)t + lc 2 (29) 



for the expression simplicity with Xq, t/o an d fi being arbitrary functions of t. 
Let us proceed to write the reduction solution p to a v = in the form 
1 , _ % „ . 1 



p= -{ C2 + 2Atanh[A(/ 2 + P)]}, A = - x / c\ + 2c Cl (30) 



with the group invariant function P = P(£, rj) and the similarity variables £, 77 being given by 

£ = + -^y 2 ftf~* + \vvotf^ -^vlft-XQ, v = yf~^-yo- (31) 

With the help of pip , the group invariant solution for other fields tp, 0, ip\, <p2 and u have the forms 

ip = # sech[A(/ 2 +P)]e~ n , (32) 

<\> = $scch[A(f 2 + P)}.r 1 / 3 e n , (33) 

^ = {-^A7^ tanh[A(/2+jP)]_ W~ l^ 04 ^/ 1 ^/ scch t A (/ 2 + P )] e ' ( 34 ) 

■tanh[A(/ 2 +P)] + ^ + ^— y ot $ + -i- sech[A(/ 2 + P)]e n , (35) 



2A/2/3 yj 18/ 12 tfUt J 2 / 3 

3//tt — 2/ 2 2 (fyot)t 1 2 M/3 1 rl/3 

18/ + 324/ 2 y + ^jm y+ 72 y ° J ~6 Xm! 

' {2PA 2 - 2c(**i + $$i)A tanh [A(/ 2 + P)] + c 2 * 2 $ 2 tanh 2 [A(/ 2 + P)] } , (36) 



2A 2 / 2 / 3 

where P = [/(£, 77), * = vj/(^ $ = 7? ) j v[( 1 = 77) and $1 = $i(£, 77) are further five group 

invariant functions while Q is related to {x, y, t} via 

n = 2f ty + 3 f p ot x £_ {lnfh |! (yot/2 /3 )t _ M__|l^ y/1 /3 _ A _ M2i. (37 ) 

Now it is straightforward to prove the following first reduction theorem. 
Theorem 2. If P is a solution of the following reduction equation, 

Ptftfif - 3P^P| - 4A 2 P«P| + 3P W P| - (1 + 4ifo)P«P* + 3P| ? = 0, (38) 

then u given by 

8A 2 P|-3P 2 -Pg-4P g P^+3P 2 c fj x 3// tt -2/ 2 2 (fy ot ) t 
U Pip/ 3 18/ 324/ 2 V 36 f 1 /^ 

+ ^ylf 4/3 - \xotf 1/3 - {P^ - APf tanh [A(/ 2 + P)]} tanh [A(/ 2 + P)] , (39) 



is a solution of the KP equation with £ and 77 being given by (|31[) , 

Proof. Substituting (|36l) into the KP equation (|T|) and vanishing the coefficients of the different powers 
of tanh[A(/ 2 + P)] = T, one can obtain seven determining equations on six group invariant functions 
U, P, VP, $, ^1 and $1. The first determining equation reads 

(2A 2 P 4 - c*$)(2A 2 P 4 + c*$) = 0, 

i.e., 

± 2A 2 P 5 - = 0, (40) 

which comes from the coefficients of T 6 . By means of (|4TJ)) . the second determining equation (the coeffi- 
cient of T 5 ) is simplified to 

c$ 2 $i +2A 2 (±*iP ? -$P ec ) =0. (41) 
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In terms of (|40|) and (|4T|) . the third determining equation becomes (by vanishing the coefficient of T 4 ) 
U = \ (i$ - ^)^V 2 - | A 2 F| + 1 (1 + 4P^)P e - 1 . (42) 



The fourth determining equation is just the equation (|38[) which is a variant form of the Boussinesq 
equation and can be read off from the coemcient of T 3 after the consideration of (|40| , ([4"T]) and (|42|) . It 
is also natural that (l38l) is just the reduction equation of the Schwarzian KP equation (fT8|) via (|30|) . 

Because the KP equation is Painleve integrable with three resonant conditions being identically 
satisfied [15], after using the relations ([40]) . (|4T|) and ([42]) . the remained three determining equations 
given by vanishing the coefficients of T 2 , T 1 and T° are all identities which have the forms 

L 2 ,i,o(P a, P c , P„, •■• , % 9„)a = (43) 

where Li, i = 2, 1, are linear operator functions of the indicated variables and 

a = left hand side of ([35]) . (44) 

Finally, the group invariant solution ([31)]) becomes ([35]) thanks to the relations (j4"0"]) . (|4ip and (14"2"]> . The 
theorem 2 is proved. 

From the expression, (|36|) (or ([39]) ). we know that when the DT related part is removed, i.e., A = 0, 
the reduction is just the known Boussinesq reduction [17j 

(U m - 6UU ( ) £ + 3U VV = (45) 

which is related to (J35J by (g2J) for A = 0. The entrance of the tanh part in (J3BJ) (i.e., ([59]) ) means the 
intrusion of an additional soliton to the Boussinesq wave. In other words, the group invariant solution 
(j39|) is an interaction solution of one soliton and one general Boussinesq wave. 

Case 2. / = 0, h ^ 0. Because the similar procedure as in Case 1, we can find the second group 
invariant solution 

P = c- 1 {c 2 + 2h 1 tanh[l(y+p')}}, 7= ^j^ , (46) 
/ // ur r 'm " (ln/i)« 3 h t gh t -2hg t 2 2F - c 2 

V, = ^Sech[ 7 ( y+P )]e , U =——y-— Sy - ^ y+——y, (47) 

= <//sech[ 7 (y + p')]e~", (48) 
V'l = |Vi- I |^[ /l *(7y+p') + 6cV'Vtanh[ 7 (y+p')]]}scch[ 7 (y+p')]e", (49) 

0i = |0i + T^[^(72/+y) + 6c^Vtanh[ 7 (y + p')]]}sech[7(y + p')]e^, (50) 

« = «' - V ~j^ {huP 1 - I2g t - hay) - ^ (VV'i + 1>W) tanh[ 7 (y + p')] 

V> ,2 </>' 2 tanh 2 [7(?/+p')] (51) 



2/i 2 7 2 

from (T u = a p = o^p = (T0 = = cr,^ = with / = 0, where u' = u'(x', t), p' = p'(x' , t), if}' = 
t), </>' = cf>'(x', t), xP[ = i/j[(x', t), ^ = ^{x', t) and 

, _ . h ty 2 . yy (M , 

x - x + l2h + T- (52) 
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Now, it is straightforward to prove the second reduction theorem of the KP equation. 
Theorem 3. If p' is a solution of the reduction equation 

2(Px'x'x'x' + Px>t)Px< ~ 2 PtP'x'P'x'x' - &l 2 Px>x>Px> - Qnh) t p'3 

+2 ((In h) tP ' - 4p' x , x , x , - Ggh- 1 ) p' x , x ,p' x , + 6p* x , ~ 6j/ a , a , = 0, (53) 



then 



Pt 4 2/2 9tP' , h tt p 12 ] g 2 h t p' - Wx'x'x' ~ 9 , 1-J& 



6p' x , 3' ex 6h 72h 2h 2 6hp' x , 2p' 2 

"-(p'htt ~ 12 flt - - 2 7 ^, x , tanh[ 7 (y + p')\ + 2 7 2 ^ 2 , tanh 2 [ 7 (j/ +p')\ (54) 



72/i 

with (|5^|) is a solution of the KP equation. 

The proof of the theorem 3 will be omitted here because it is similar to the proof of the theorem 2. 

From the expression, (f31)|) (or (|3"9"10 . we know that when the DT related part is removed (c = 0), the 
reduction will be reduced back the known KdV type reduction |17j 

(Uix>x>x> - GUiUtxA, + U lxH + ^-U lx , + ^- = (55) 
v x In 12n 



with 



Ul = U -^2h- + -6h-2h?- (56) 



The appearance of the tanh part in (|54p implies that the group invariant solution is an interaction solution 
of one soliton and one general KdV type wave. 

Case 3. / = h = 0. In this case, the group invariant conditions a u — cr p = ay, = — = a c / >1 = 
with / = h = possess the solution 



p = - {c 2 - 2. 9 /3 tanh\J3(x + q)]} , /3 = ^ +2cCl ; (57) 
c 2g 

ij = rsech[l3(x + q)]e Ul , ui = ~(V9t + - 3c 2 ), (58) 
= s sech[/3(a; + g)]e~ Wl , (59) 

(CT^ S \ 
— tanh[/3(x + 9 )] + n 1 sech[/3(x + q)]e" 1 , (60) 

— ■ tanh[/?(x + <?)] + Si sech[/3(a; + q)^ 1 , (61) 
25P J 

2 2 2 

u = ui + +rsi)tanh[/3(a; + q)] + tanh 2 \ft(x + q)] - f^( x + l) (62) 

with the group invariant functions q = q(y, t), r = r(y, i), s = s(y, t), n = r\{y, t), s\ = Si(y, t) and 
ui = m(y, t). 

Similar to the last two cases, substituting (|6"2"|) into the KP equation leads to the following third 
reduction theorem. 

Theorem 4. If q is a solution of the reduction equation 

q vv = 0, (63) 



11 



then 

u = ~<?t + \ql - \p - f g {x + q) + 2/3 2 tanh 2 [/3(.T + q)] (64) 

is a solution of the KP equation. 

It is clear that the reduction solution (|64|) is corresponding to the usual single soliton solution (for 
g = constant and q being a linear function of y and t) or the single solitary wave with nonhomogeneous 
background for g being not a constant. 



5 Interactions among different nonlinear waves 

From the reduction theorems 2, 3 and 4 of the last section, one can find some important aspects: (1). 
The DT related symmetries can be used to find new symmetry reductions which are the generalizations 
and will be reduced back to those of without DT transformations. (2). The essential role of the DT is 
to add one soliton no matter what the seed solution is. (3). The interactions between solitons and other 
nonlinear waves can be analytically found. The theorem 2 and theorem 3 display the interactions among 
solitary waves and Boussinesq waves and KdV waves respectively while the third reduction is related to 
the single solitary waves with homogeneous or nonhomogeneous background. 

For concreteness, we write down two further explicit solutions of the first two reduction equations. 

For the first reduction equation (|38[) . its cnoidal wave solution is given by, 

R 2 X1 -A 2 R i -R + C 1 R 2 -C 2 R 3 = 7 R = P Xl , P(£, r?) = P(£ + ar]) = Pfa), (65) 

with arbitrary constants a, C\ and C%. The solution of (|65|) can be explicitly expressed by Jacobi elliptic 
functions, say, 



R = ^7 — 52' S = sn (^V / ( r 'i -r 3 )r 2 Axi,m) , m= J'^p- — (66) 



where S is the Jacobi elliptic function with the modulus m and the arbitrary constants {Ci, C%, A} 
have been re-expressed by 

Ci = r^+nrs+rara ^ = n + r 2 + r 3 ^ = 1_ 

nr 2 r 3 nr 2 r 3 r x r 2 r 3 



Correspondingly, the field P has the form 
2r 3 



P = 



E * 1 S ^ n -r 3 + 2r lS ^ W^- 3 m ) ~ EF N n - r's + 1^ ™ 



(68) 



&\/r 2 {r x - r 3 ) 

where Ep and E v are the first and third incomplete elliptic integrals. It is clear that the exact solution 
(|39[) denotes the interaction between the Boussinesq-type cnoidal periodic wave and soliton by taking 
the arbitrary functions /, y and Xq as / = l/i>i, Xq — v 2 t, y = ^3(^ + ^0) with arbitrary constants 
vi, v 2 , v 3 and t . 

For the second reduction equation (|53[) given in theorem 3, we consider only h = 1, g = Co = constant 
case. For the constant h and g case, its cnoidal wave solution is given by, 

r 2 yi - 4 7 V - 1 - 6c R - CV - C 2 r 3 =0, r = p' yi , p'{x', t) = r{x' + at) = r(yi), (69) 
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where C\, C 2 and a are arbitrary constants. 

The general solution of (JSHJ) can be written as 



(n -n)nSl -r 4 (n -r 3 ) _ / V( r i - r 3 (r 2 - r 4 ) \ _ /(n -r 4 )(r 2 -r 3 ) . . 

r = Ti w2 — > Si = sn y u m , m = W- r (70) 

^3 ^ n + [n - Ujbf \ 2rir 2 r 3 r 4 y y (r 2 - r 4 )(r x - r 3 ) 

if the arbitrary constants 7, Co, Ci and C 2 is re- written as 

^ _ n + r 2 + r 3 + r 4 _ r x (r 2 + r 3 + r 4 ) + r 2 (r 3 + r 4 ) + r 3 r 4 

Oi — , 2 — , 

?"ir 2 r 3 r 4 rir 2 r 3 r 4 

, 1 1/1 1 1 1 \ 

1 2 = - A , c = -- _ + _ + _ + _ . (71) 

4rir 2 r 3 r 4 6 \r x r 2 r 3 r 4 / 



Consequently, the field p' has the form 



V 



4rir 2 r 3 r 4 



(r 3 - r 4 )^ ( Si, — — ,m I - r 3 ^ F (5,m) 

ri - r 3 



(72) 



(r 2 - r 4 )(n - r 3 ) 
6 Discussion and Summary 

To sum up, the non-local symmetry of KP equation related to the DT can be localized if five 
potentials, the spectral function ip, the adjoint spectral function tp, the x-derivatives of the spectral 
functions = tpxi ipi = ^xt } and the singularity manifold function p = J -0</>dx, are introduced. For 
the enlarged system {u, ip, cj>, <f>i, p} the generalized Lie point symmetries includes the DT as its 
special ones. Using the Lie's first theorem to the Lie point symmetries of the enlarged system, one can 
find that the first type of DT and the second kind of DT possess the same infinitesimal form but with 
different group parameters. From the infinitesimal form of the DT, is is found that the DT of the KP 
equation is nothing but the Mobious transformation of its Schwarzian form. This conclusion is true for 
the KdV equation [TU] and may also correct for other integrable systems [T5] . 

Usually, to find the new solutions via DT, one has to solve the spectral problem first. However, for 
non-constant and non-soliton potential (seed solution) it is very difficult to solve the spectral problem. 
Fortunately, In this paper, it is found that it is not necessary to solve the spectral problem directly. The 
simple alternative way is to find the symmetry reduction method with the DT related symmetries. The 
result shows that if we take the Boussinesq type and the KdV type waves (the general solutions of the 
Boussinesq and KdV reductions of the KP equation) as seed solutions, the essential and unique role of 
the DT is to offer an additional soliton. 

The KP equation was widely used in various physics fields such as the surface and internal oceanic 
waves, nolinear nonlinear optics, plasma physics, ferromagnetics, Bose-Einstein condensation, string the- 
ory and so on. Though many kind of solutions of the KP equation have been found by many authors via 
different methods, however, it is difficult interaction solutions among different kinds of nonlinear waves. 
Especially, the soliton solutions can be used to described various nonlinear phenomena. For instance, 
the KP's line soliton solutions can be used to describe the propagating tsunami wave, on the other 
hand, there are many other oceanic surface waves such as the periodic waves which may be described 
by the Jacobi elliptic periodic solutions (cnoidal waves) of the KP equation. Then how can we find the 
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interaction solutions among the line solitons and the cnoidal waves? In this paper, two special explicit 
such kinds of solutions are obtained as shown in (|3T)1) with (1551) and with ([72)1 respectively. 

The following topics will be discussed in the future series research works: (1). The method used here 
can be developed and/or extended to other kinds of nonlocal symmetries such as those of obtained from 
Backlund transformation [TO], the bilinear forms and negative hierarchies [191 120) . the nonlinearizations 
[21] . self-consistent sources [22] and the Painleve analysis [15.. (2). It is shown here that to find the 
classical symmetry reductions of the enlarged system is equivalent to get the DTs with the seed solutions 
being taken as the classical symmetry reductions of the original nonlinear system. The method is valid and 
the conclusion is correct for other kind of complicated seed solutions such as the nonclassical symmetry 
reductions. (3). The method should and can be applied to other kind of integrable systems, especially for 
supersymmetric models and discrete ones, to find interaction solutions among different kinds of nonlinear 
waves. (4). Other types of well known effective methods can also be extended to study the interactions 
among different kinds of nonlinear excitations. 

The work was sponsored by the National Natural Science Foundation of China (No. 11175092), 
Scientific Research Fund of Zhejiang Provincial Education Department under Grant No. Y201017148, 
and K. C. Wong Magna Fund in Ningbo University. 
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